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We analyze a randomly perturbed quantum version of the baker's transforma- 
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£f~^ ' tion, a prototype of an area-conserving chaotic map. By numerically simulating 

ON ■ the perturbed evolution, we estimate the information needed to follow a perturbed 

Hilbert-space vector in time. We find that the Landauer erasure cost associated with 
this information grows very rapidly and becomes much larger than the maximum 
statistical entropy given by the logarithm of the dimension of Hilbert space. The 
quantum baker's map thus displays a hypersensitivity to perturbations that is analo- 
gous to behavior found earlier in the classical case. This hypersensitivity characterizes 
O ■ "quantum chaos" in a way that is directly relevant to statistical physics. 



Great progress has been made in studying manifestations of chaos in quantum systems 
yet there still is controversy as to whether quantum chaos exists at all || [3[ . A chief reason 
for this is that the most important characteristic of classical chaotic systems — exponential 
divergence of trajectories starting at arbitrarily close initial points in phase space — is absent 
from quantum systems simply because the existence of a quantum scale makes meaningless 
the concept of two arbitrarily close points in phase space. 

Any attempt to find exponential divergence of trajectories of Hilbert-space vectors 
founders immediately, because the linear Schrodinger equation, with its unitary evolution, 
preserves Hilbert-space inner products. Yet the unitary linear evolution of the Schrodinger 
equation must be irrelevant to the issue of quantum chaos 0, since any Hamiltonian 
classical chaotic system can be described by an analogous area-conserving linear Liouville 



"Submitted to Physical Review Letters 



1 



equation. Two probability distributions in classical phase space — we call these phase-space 
patterns — that are initially close together (in terms of an overlap integral) stay close to- 
gether forever, because the Liouville equation is area conserving. This leads to the question 
we address in the present article: Are there manifestations of Hamiltonian classical chaos 
in the Liouville equation that are also present in quantum mechanics? 

In an earlier paper |4[] we analyzed a prototype of an area-conserving chaotic map, 
the baker's transformation ||, in the Liouville representation; i. e., we focused on phase- 
space patterns instead of on single trajectories. Our analysis was guided by the question 
H of how available work decreases with time when the baker's map is subjected to area- 
conserving random perturbations. An area-conserving abstract mapping corresponds to 
an energy- conserving phase-space system, so we identify two negative contributions to 
free energy. The conventional one is ordinary entropy, which measures how incomplete 
knowledge about a system reduces our ability to extract work. The other contribution 
arises from Landauer's principle |7|, [J that there is an unavoidable energy cost of /csTln2 
connected with the erasure of one bit of information. It follows from Landauer's principle 
that the information, quantified by algorithmic information [[|, needed to give a complete 
description of a system state also reduces the amount of available work and thus should 



make a further negative contribution to free energy 0, [TT[ . 

In our earlier paper M we compared two strategies for preserving the ability to extract 
work from a system. The first strategy is to keep track of the perturbed phase-space pattern 
in fine-grained detail, in an attempt to preserve the work inherent in the initial condition. 
The second strategy, which we call coarse graining, is to average over the perturbation and 
to put up with the resulting ordinary entropy increase. We found, for the perturbed baker's 
map, that the information needed to implement the first strategy is overwhelmingly larger 
than the entropy increase of the second strategy. This means that the free-energy cost 
of tracking the perturbed pattern in fine-grained detail is enormous and far greater than 
the cost of the entropy increase that results from coarse graining. We conjecture that this 
hypersensitivity to perturbations is a general feature of perturbed classical chaotic systems, 
and we regard it as the desired manifestation of classical chaos in the Liouville equation. 

In the present paper we compare the two strategies for preserving work in the case of a 



quantum system, a quantum version of the baker's map [fL2jl . Using numerical simulation, 
we find essentially the same behavior as in the classical case, as was suggested using heuristic 
arguments in Ref. ||. The quantum baker's map displays hypersensitivity to perturbations 
and thus can be said to exhibit quantum chaos. 

The concept of algorithmic information has been used before to investigate quantum 
chaos H |13| . If one defines a chaotic system as one where the algorithmic information 
needed to predict a single (unperturbed) trajectory grows linearly with time (or number 
of steps) PH , then there is classical chaos, but no quantum chaos @. Our approach, by 
focusing on patterns in phase space instead of trajectories, uses a framework where classical 
and quantum mechanics can be treated on analogous footings. Moreover, since Landauer's 
principle gives information an explicit physical meaning by connecting it to available work, 
our characterization in terms of hypersensitivity to perturbations is directly relevant to 
statistical physics. 

The classical baker's transformation maps the unit square < q,p < 1 onto itself 
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according to 

f:{q,p)^(2q-[2q\,(p+[2q\)/2), (1) 

where the square brackets denote the integer part. There is no unique way to quantize a 
classical map. Here we adopt a quantized baker's map introduced by Balazs and Voros 



12| and put in more symmetrical form by Saraceno |15]. Position and momentum space 
are discretized, placing the lattice sites at half-integer values qj = pj = (j + 1/2) /2N 
for j = 0, . . . , 2N — 1. The dimension 2N of Hilbert space is assumed to be even. For 
consistency of units, let the quantum scale on phase space be 2irh = 1/2N. Position and 
momentum basis kets are denoted by \qj) and \pj). A transformation between these two 
bases is performed by the operator GW, defined by the matrix elements 



(G 2N ) jk = {pj\q k ) = e- ip ^ n . (2) 

The quantum baker's map is now defined by the matrix 

B — &2N { n n J ' (^) 



where, as throughout this article, matrix elements and vector coordinates are given relative 
to the position basis. 

The perturbation operator we use is constructed to resemble the type of perturbation 
used in our previous work [|] on the classical baker's map. We partition phase space into 
an even number 2N C = 2N/w c of congruent perturbation cells, where 2N C and w c > 2 are 
integral divisors of 2N. In the following we use perturbation cells that are vertical stripes 
extending over the entire p range. Then each perturbation cell contains w c g-eigenstates. 
A perturbation operator that perturbs each perturbation cell independently has the form 
of a matrix with zero elements everywhere except for 2N C square blocks of size w c on the 
diagonal. We choose these square blocks so that they correspond to a shift in the p direction 
in a w c - dimensional subspace. Let the momentum shift in the nth perturbation cell (n = 
0, . . . , 2N C — 1) be ai n , where the real number a is the magnitude of the momentum shift 
and i n G {—1, 1}. The symmetry condition i2N c -n-i — in (n = 0, . . . , N c — 1) avoids rapid 
oscillations and thus ensures similarity to the classical case. The perturbation operator 
;io,...,ijv c _i is defined by 

(Ua;io,— ,iisi c -\)jk $jk& ^ k j (4) 

where O — 0, <pk — 4>k-i + a ^n(fc) {k = 1, . . . , 2N — 1) with n(k) = [k/w c ]. The parameter 
a characterizes the "strength" of the perturbation, whereas w c /2N = 1/2N C is the area of 
the perturbation cells. 

A perturbed time step consists of first applying the unperturbed time-evolution operator 
B, followed by a perturbation operator U a;i0; _ t i N _ x with i , . . . ,iN c ~i £ { — 1> 1} chosen at 
random, a being fixed. We thus allow for a different perturbation at each step, in contrast 
to Ref. [0 where a particular perturbed evolution operator was applied repeatedly. After 
n time steps, the number of different perturbation sequences — or histories — is 2 nNc . 

Our specialization to vertically striped perturbation cells involves no restriction relative 
to our work on the classical baker's map. There we allowed for (2 m w c /2N) x 2~ m rectangular 
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perturbation cells, which are the image of vertical stripes under m applications of the baker's 
map, 2 m < 2N/w c = 2N C . Likewise, the perturbation operator for "rectangular" quantum- 
mechanical perturbation cells is U' = B m UB~ m where U is given by Eq. (^). Using U with 
initial state \ipo) is equivalent to using U' with initial state B m \ip ). The freedom to choose 
m — we use m = for vertical stripes — is the same as the classical freedom to choose the 
initial position of the "decimal point" in the symbolic representation of the baker's map. 

As a preliminary step, we show that perturbed evolution leads after several steps to 
an ensemble of vectors that is similar to an ensemble of vectors distributed randomly on 
Hilbert space. A useful criterion for determining the randomness of an ensemble of vectors 

is based on the moments of the quantity Wt 



2N-X 

W(m = - E |<^>| 2 log(|(g#>| 2 ), (5) 

which we call PU-entropy to distinguish it from ordinary entropy (log denotes the binary 
logarithm, as throughout this paper). For random vectors in 2N = 16-dimensional Hilbert 
space, the mean and standard deviation of the IV-entropy are given by W = 3.434 ± 0.178 
bits, a result obtained on a computer by calculating H^(| , 0)) for a large number of vectors 
chosen at random from an ensemble distributed uniformly over Hilbert space fTT 



(This mean value agrees with the exact formula for the mean value W = \${2N + 
1) — ^(2)]/ In 2, where \I/ is the digamma function.) We compare the result for random 
vectors with the first two moments of the jy-entropy for an ensemble of vectors created 
by the perturbed baker's map in a 16-dimensional Hilbert space. Choosing 2N C = 2 and 
a = 0.025 = 0.4/2iV, we create an ensemble of approximately 20 000 perturbed vectors 
by applying different randomly chosen histories for n = 15 perturbed time steps to the 
initial vector \ipo) = \ps). We find for the first two moments of the jy-entropy the values 
W = 3.428 d= 0.184 bits, very close to the moments for a random sample. 

As a second check of randomness, we calculate the distribution of Hilbert-space an- 
gles 9 = cos - 1 ( | (-i/^' | ) between vectors \ip) and \ip') that have evolved under the same 
perturbed quantum baker's map applied to the same initial state as in the previous ex- 
ample. We compute the Hilbert-space angle between each pair of vectors in each of three 
ensembles of approximately 16 000 vectors created by applying different randomly chosen 
perturbation histories for 15, 23, and 31 steps. In addition, we compute the Hilbert-space 
angle between each pair of the 2 10 vectors after 10 steps. The resulting distributions of 
Hilbert-space angles are displayed in Fig. |l|. After 31 steps the closest pair of vectors is 
27.1° apart, a striking demonstration of the "size" of 16-dimensional Hilbert space, i. e., of 
how many widely separated vectors Hilbert space can accommodate even for a relatively 
small dimension. For comparison, Fig. |l]also shows the distribution f(6) = 30(sin6') 29 cos 9 
of Hilbert-space angles for a set of random vectors. These results show clearly how the 
ensemble is randomized by the perturbed quantum baker's map. 

We proceed now to compare the two strategies for extracting work outlined above — 
coarse graining versus following the evolved vector in fine-grained detail. We estimate 
the conditional algorithmic information AI needed — given background information — to 
specify a typical perturbed vector after n steps and compare it to the increase in ordinary 
entropy AH that results from averaging over the perturbation. Our first example uses, as 
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Figure 1: Distribution of Hilbert-space angles for vectors evolving under the perturbed 
quantum baker's map, shown for different numbers of perturbed time steps n. For com- 
parison, the distribution for random vectors is also shown. The dimension of Hilbert space 
is 2N = 16, the number of perturbation cells is 2N C = 2, the perturbation strength is 
a = 0.025 = 0.4/2 A/", and the initial vector is |Vo) = \ps)- 



before, a 2N = 16-dimensional Hilbert space, partitioned into 2N C = 2 vertically striped 
perturbation cells. We choose a fixed perturbation amplitude a = 0.025 = 0A/2N and 
an initial pure state |^ ) — \ps), i- e -> a momentum eigenstate, which corresponds to a 
horizontal stripe in the unit square. This perturbation can be described completely by 
giving one bit per step, to specify which of the two possible perturbation operators U a -i 
and U a --i is applied. If the logarithmic term |IT| that keeps track of the number of steps 
n is neglected, this sets an upper bound on the information AI. This upper bound is 
realized only if two different histories of perturbed time steps always lead to two different 
vectors at some level of resolution on Hilbert space. We choose a resolution that regards 
two vectors as different if their Hilbert-space angle exceeds 89 = tt/50 = 3.6° (| (ip'\ip) | 
smaller than 0.998). By comparing numerically all possible histories, we find that, through 
15 perturbed time steps, all trajectories lead to distinguishable vectors. Figure |2| shows 
the resulting linear increase in the information AI. 

Figure |2| also shows the ordinary entropy increase AH, obtained by determining the 
entropy of the density matrix that results from averaging over all possible histories. It 
can be seen that AI is always larger than AH. Indeed, AH saturates at the value 
AH mSuX = log 2N = 4 bits, the logarithm of the dimension of Hilbert space, whereas AI is 
only limited by A/ max ~ — 2(2iV — 1) \og(59/2) ~ 150 bits, which is the logarithm of the 
number of different vectors Hilbert space can accommodate || [HJ . Whereas Aif max grows 
logarithmically with the dimension 2N of Hilbert space, the maximum information A/ max 
grows linearly with 2N and is enormous for macroscopic systems. 

As in the classical case [|j] , the information AI grows more dramatically when the num- 
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Figure 2: Conditional algorithmic information AI needed to track a vector evolving under 
the perturbed quantum baker's map, compared to the increase in ordinary entropy AH 
that results from averaging over the perturbation. The parameters are the same as in 
Fig. 1. 
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Figure 3: AI and AH as in Fig. 2, but with parameters 2N = 64, 2N C = 16, a = 0.05 ~ 
3/2iV, and |^ ) = B~ 5 \ Pl ). 
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ber of perturbation cells is large. Figure |3] displays results for a 64-dimensional Hilbert space 
with 16 vertically striped perturbation cells, each containing four position eigenstates. The 
perturbation strength is a = 0.05 ~ 3/2N, and the initial state is l^o) — B~ 5 \pi), a state 
whose image under B has negligible support outside the leftmost perturbation cell. This 
means that, in order to describe the perturbed state after the first time step, in the pertur- 
bation operator U a . io> ___ ji7 only the sign i referring to the leftmost perturbation cell must be 
specified. Since B 2 \ip ), B 3 \ip ), and B^ipo) extend over 2, 4, and 8 perturbation cells, we 
expect the number of bits needed to specify the perturbed state to grow as X^=o ^ = 2 n — 1 
until the state extends over all perturbation cells. This behavior is verified in Fig. |3] using 
the same method as for Fig. ^|. 

Given these results and those of our previous paper J3J, we have demonstrated similar 
hypersensitivity to perturbation in a classically chaotic system and its quantum analogue. 
In both cases the large information needed to track the perturbed evolution is due to the 
large number of possible ways to perturb a state || |21J . In the classical domain, chaos opens 
up the large space of possibilities — phase-space patterns with structure on finer and finer 
scales. Quantum mechanics operates inherently in an enormous space of possibilities — the 
pure states on Hilbert space. Hypersensitivity to perturbations means that more work can 
be extracted by coarse graining than by following the perturbed evolution in fine-grained 
detail. Our results provide a motivation for coarse graining and thus an explanation of the 
second law of thermodynamics. 

RS acknowledges the support of a fellowship from the Deutsche Forschungsgemeinschaft. 
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